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Abstract 

Low energy QCD (below 2 GeV) is a region of resonance dynamics, sometimes 
lacking satisfactory description as compared to precision of available experimental 
data. Hadronic r decays offer a probe for such energy regime. In general, predictions 
for decays are model dependent, with parameters fitted to experimental results. Pa- 
rameterizations differ by amount of assumptions and theoretical requirements taken 
into account. Both model distributions and acquired data samples used for fits are 
results of complex effort. 

In this paper, we investigate main parameterizations of r decay matrix elements 
for the one- and three-prong channels of three-pion r decays. Differences in analyt¬ 
ical forms of the currents and resulting distributions used for comparison with the 
experimental data are studied. We use invariant mass spectra of all possible pion 
pairs and the whole three-pion system. Also three-dimensional histograms spanned 
over all distinct squared invariant masses are used to represent results of models and 
experimental data. 

We present distributions from TAUDLA Monte Carlo generation and semi-analytical 
calculation. These are necessary steps in development for fitting in as model-independent 
way as possible, and to explore multi-million event experimental data samples. This 
includes response of distributions to model variants, and/or numerical values of pa¬ 
rameters. Interference effects of currents parts are also studied. 

For technical purposes, weighted events are introduced. Even though we focus on 
Stti/t modes, technical aspects of our study are relevant for all r decay modes into 
three hadrons. 
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1 Introduction 


One aspect of phenomenological work is to represent theory and experimental efforts in 
form of Monte Carlo generators. Then, user should be able to obtain any distribution 
of interest. All details of the detector response can be taken into account if needed. 
Introduction to the tasks of modeling with the help of Monte Carlo of low energy hadronic 
interactions can be found, e.g., in [1], and will not be repeated here. We will concentrate 
on T lepton decays. 

Let us point out that there are various aspects of the work for construction of Monte 
Carlo, which need to be combined. On the theoretical side, assumptions have to be trans¬ 
lated into distributions, which can be confronted with the data. This task is easy, if un¬ 
derlying theory is well-established and its calculation methods based, e.g., on perturbation 
expansions, are converging sufficiently fast. For theorists, it would be ideal if experimental 
results were represented in a form of background-subtracted and detector corrected distri¬ 
butions, sufficient to constrain all constants (masses, widths, couplings) introduced by the 
theory (or model) used as a basis for the calculation. 

In practice, such a task is often far from being straightforward. Assumptions behind 
theoretical calculations are often not well-established and the number of phenomenological 
constants introduced may be too large to be constrained by experimental data. 

From the point of view of Monte Carlo techniques, generating processes such as r 
decays into Svri/ is rather simple. However, necessary for that purpose, hadronic currents 
span distributions over eight-dimensional space and are the result of a massive effort. 
Taking into account Lorentz invariance and properties of weak couplings of r-lepton to 
intermediate virtual W boson, predictions still require four complex scalar functions F) of 
three variables each (see the dehnition later in the text). Such parameterization of matrix 
elements is commonly used by r decays Monte Carlo generators, e.g., by TAUOLA; and since 
its beginning Pi. 

In principle, all properties of the matrix element can be constrained from experimental 
data. All F, can be htted in a model independent way as proposed in [1]. In practice, 
this is highly non-trivial. Even if necessary for that purpose, direction of r neutrino 
can be reconstructed, one has to measure from the data at least 7 distributions over three 
dimensions with sufficient detail. So far, this was never achieved without partial integration 
over some of the phase-space dimensions. This practical aspect of the phenomenology work 
should be separated from the previous two (theoretical and experimental ones), and at the 
same time offer comfortable and flexible usage. 

In this paper, we investigate differences and physics motivations for the presently used 
in TAUOLA hadronic currents of the r —)■ Sttz/,- decay channels. Let us recall from Ref. [5] 
comparisons of different parameterizations for case of the decay mode r —>■ 27ri>r. In that 
case, all results are encapsulated in a single one-dimensional distribution. Our analysis, as 
we will see, must rely at least on three-dimensional distribution. We compare its efficiency 
to constrain model parameters hrst with the method when one-dimensional histograms are 
used only. To evaluate full size of the model differences, integral of —1| = | —1| 

over event samples will be used. M^odeh MmodeV denote matrix elements for the two 
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compared model variants. For completeness, we study the case when only total width of 
the channel is used to constrain the model. 

Theoretical principles embedded in the models and remaining freedom for introducing 
new couplings and/or intermediate resonances has to be studied at the same time. We 
have to keep in mind that additional contributions to the currents can modify shapes of the 
distributions in an nonphysical/uncontrolled way and in unexpected regions of the phase- 
space. This is especially important when dimensionality of experimental distributions is 
smaller than that of the decay phase-space. 

The paper is organized as follows. In Section [2l we present the list of currents used 
for evaluation of distributions. Detailed description is delegated to Sections El 01 [5l and 
ini analytic forms of the currents are given. In Section [71 we present basic features and 
motivations behind models used in data htting. Section |8] consists of introduction to 
numerical comparisons and commonly used representations of experimental data employed 
in fits. Section El is devoted to models used by the CLEO collaboration. In the subsections, 
numerical results and comparisons are presented. Section HU] describes further variants of 
currents that were never used outside the experimental collaboration but their existence is 
mentioned in Ref. |6]. Two main CLEO models are then compared with the BaBar model in 
Section [TTl Section [T2]provides comparisons between default TAUOLA model and two models 
based on the Resonance Chiral Lagrangian (RCliL) [7]. For this task, we use mainly one¬ 
dimensional histograms; results of measurements in such a form are available not only from 
the CLEO but from the BaBar as well. Such distributions were used to constrain model 
parameters in many of our studies. Only CLEO used the three-dimensional distributions 
as an input for fits. 

In Section [13] contributions from different parts of currents and their interferences are 
investigated. We discuss possible consequences of using model assumptions and limited 
data samples in Section [TT] Summary, Section [TH] closes the paper. In Appendix [^ we 
present arrangements for semi-analytical calculations which can be used for hts into three- 
dimensional histograms. In Appendix [B] numerical comparison of results from different 
models which are of smaller importance are given. Lengthy histograms of three-dimensional 
form, are given in Appendix O 

2 List of currents 

In TAUOLA Monte Carlo generator mu multiple versions of hadronic currents for the two 
channels of r —)■ Sttz/,- decay are available. They are the result of theorists effort but also of 
extensive work of the experimental collaborations. Over time, currents became available 
for use outside of the experiments. Unfortunately, sometimes these currents are poorly 
documented. We plan to explicitly present differences in formulae and values of numerical 
constants for the main currents available for TAUOLA. For others, we will skip some details 
as they are of lesser importance or are documented in quoted publications. 

The following variants of the currents will be discussed: 

1. CLEO publ. - current described in ref. [6]. This parameterization was developed 
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for data of the r Ur decays. Details will be discussed in Sections H] and [H 

2. TAUOLA CLEO - current implemented in TAUOLA. Code of this current was ob¬ 
tained from the CLEO collaboration and distributed with TAUOLA (all its versions 
since Ref. m)- See Sections 0] and |9] for details. 

3. TAUOLA CLEO isospin intricate - this is a variant where isospin rotation from 
the 7r°7r°7r“ to the 7r“7r“7r’'' channel was performed for intermediate resonances. Nu¬ 
merical constants remained unmodihed. Code was obtained from the CLEO collab¬ 
oration and distributed with TAUOLA (all versions since [8]) but it was never activ^. 
See Sections 0] and [9] for details. 

4. TAUOLA BaBar - parameterization numerically equivalent to the one used in the 
BaBar collaboration for basic simulations. Details are given in Sections |5] and [11] 

5. TAUOLA RChL 2012 - version of RChL model introduced in TAUOLA in 2012. This 
model was based on theoretical consideration of that time and hts to invariant mass 
distributions of 7r“7r“7r+ and systems [13]. For details see Sections 131 and [T^ 

6. TAUOLA RChL - this is the model motivated by further theoretical consideration, 

and comparison with experimental data also for invariant mass, Ref. |14j . 

7. TAUOLA CPC - outdated parameterization of Monte Carlo generator. Founda¬ 
tion of many technical benchmarks for phase-space and other verihcations which are 
documented in [3]. Also, test distributions documented in Ref. are available. 

8. Pythia CLEO - the parameterization of Ref. |6] as implemented in Pythia [15]. Cur¬ 
rent nearly identical to the TAUOLA CLEO isospin intricate is used. Currents 
of Pythia differ only by numerical constants. In the present paper, we use numerical 
results of Pythia version 8.201. This current will be discussed in Appendix iBl 

At present, TAUOLA CLEO is the default choice which should give the same re¬ 
sult regardless of which version of TAUOLA is used. It remained unmodihed for at least 14 
years. It is identical for FORTRAN and C++ implementations. Differences in distributions 
between CLEO pnbL, TAUOLA CLEO, TAUOLA CLEO isospin intricate and 
Pythia CLEO are presented with all technical and numerical aspects. For completeness, 
we present numerical results from TAUOLA RChL. This last current is aimed to rep¬ 
resent data with good precision, and at the same time improve relation with present day 
theoretical calculations. Numerical results for other variants of the currents are given in 
less detail or are not presented at all. 

^ This parameterization is missing validation with experimental publication; only conference contribu¬ 
tions [iQlITT] announce its existence, while PhD thesis [12] fully describes even more elaborated variant. 
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3 Common mathematical functions 


All currents parameterizations use the same basic formulae. Equation ([I]) define^ hadronic 
current. It is common for all three-scalars decay channels of TAUOLA and with its help 
matrix element of r decays can be calculated in a straightforward manner Pi, 


JM ^TI^^Ci{p2-p3)''Fi+C2{p3-pi)''F2+C3{pi-p2)‘'F3] 


+C4q^F4 - 




P 1 P 2 P 3 F 5 }, 


( 1 ) 


where — Q^^Qv/Q^ denotes the transverse projector, and = (pi -|-p 2 +^ 3 )^ 

is the momentum of the hadronic systemic. The decay products are ordered 7 r“ 7 r“ 7 r’'' for the 
three-prong channel and for the one-prong channel and the pions four-momentc 0 

are denoted as pi, p 2 and pa, respectively. The same ordering is used in further equations 
also for masses (ruj). The yp„ is the Levi-Civita symbol. In equations of this and the 
following sections we use notation: s, = (pj + Pk)"^ where i 7 ^ j 7 ^ k 7 ^ i. Constants: ci, 
C 2 , etc. are Clebsch-Gordan coefficients, dehned specihcally for particular hadronic current 
used. 

Eqs. (E]), ([3]) and (H]) describe Breit-Wigner functions that are later used in dehnition 
of form factors: 


P{S,mi,m2) 


^/{S - {mi + m 2 ) 2 )(S' - (mi - 1712)“^) 

Vs 


Tl —wave {S, M, r, mi,1712, L) 


/ P{S,mi,m2) 


BW{S, M,r,mi,m 2 ,L) 


S - M 2 - iMTL^^aveiS, M, r, mi, m 2 , L)' 


( 2 ) 

(3) 

(4) 


These are typical building blocks useful for hadronic currents parameterizations, such 
as of Gounaris-Sakurai parametrization [16] for p —)■ tttt. 


4 Hadronic current in models of the CLEO category 

Let us describe the default, hadronic currents of the CLEO modeling - TAUOLA CLEO. 
In eq. (El), ® the analytic form of the hadronic current for the r —)■ decay channel 

is presented^ Its functional form is extracted from TAUOLA cod^ but has a form exactly 

^ We use five form factors instead of four imposed by Lorenz invariance for practical purpose. In principle 
can be represented as linear combination of contributions to Fi and F 2 . 

^In TAUOLA code defined as PAA(4). 

^ In TAUOLA code represented as PIMl (4), PIM2(4), PIM3(4). 

® The vp^PiP 2 P 3 is coded in subroutine PROD5. 

® Constants: f3,l3i,fi2 etc. are coded in TAUOLA as BET, BTl, BT2, etc. 

^Such approach checks if over the years any modifications were (intentionally or not) introduced. 
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as described in appendix A of Ref. [6]. Notations differ slightly. This form of hadronic 
current is also used for the r —)■ 27T~7r~^i>T- decay channel in default TAUOLA setup. It will 
be discussed later why such a choice is justihed. In PDG tables Ref. ra , some mesons 
have the same name but different mass. In our paper; ai(1260), a (or /o(500)), /2(1270), 
/o(1370), p(770), p'(1450), K*{892) are used. Finally, meson a']^(1700) as in Ref. [6]. Masses 
in names will be skipped from this point on. 


G = 


Ml 


Ml 


- Ml - ■ WGA{Q^) q2 _ • WGA{Q^) 

/ 3 i ■ BW{si,Mp,Tp,m 2 ,m 3 ,l) + P2 ■ Mp-, Pp/, m2, m3,1) 

(s3 - m§) - (si - ml) 


-da 


• BW{s2,Mp,rp,m3,mi,l) 


g -.jp, i p, 

n2^ _ I _ ^2\ 


{S 3 - ml) - (si - ml) 


■ ds 


3 “ 

iQ^ + S3 - rn|)( 2 m§ + 2ml - S3) 
2 


BW{s2,Mp>,Tp>,m3,mi, 1) 

BW{s3,Mf^,Tf^,mi,m2,2) 


+ d6 • 2 ■ BtP(s3,Mo.,rcr,mi,m2,0) 

+ d7 • ^ • BW{s 3, Mfg,rf^,mi,m2,0)'^ , 


(5) 


F 2 has the same functional form as Fi. The only difference is interchange for its 
arguments indices 1 and 2 in eq. CD. and that constant C 2 has opposite sign to ci. Note, 
that (3 is always set to 0. We suspect that it was introduced by the CLEO collaboration 
for studies of influence. 


F 3 = 


Ml 

ai 


Ml 


- Ml - ■ WGAiQ^) Q2 _ Ml - ■ WGd(Q2) 

da • ~ 3 ~ • i?iT(si,Mp,Pp,m2,m3,1) 

^ . bvP(s 2, Mp, Pp, m3, mi, R 

d 4 - ~ "^ 2 ) - (sa - mi) .B^^^^^Mp,,Tp,,m 2 ,m 3 A) 


(s3 - ml) - (si - ml) 


■ BtP(s2, APp-,Pp/,m3,mi, 1) 


- da • ^ ■ BW{s3,Mf^,rf^,mi,m2,2) ). 


( 6 ) 


5 


















Form factors F 4 and F 5 are set to zerc@. Note, that normalization factor used in oi and 
a'l width is the same. The dependence of the Oi width in its propagator is given by 
the formulate ([7]). It was obtained from analysis performed by the CLEO collaboration [ 6 ]. 
Contributions from the three main oi decay channels are taken into account. Complicated 
form is determined by the decay channels: as ai virtuality gets larger thresholds are crossed, 
allowing more decay channels to open, therefore changing the dependence of the effective 
width. is given in GeV^ units. In the equation below: = 0.2384^ and Ck* = 

4 . 76212 ^ 3 .. 


WGA{Q^) = 


1 " 

■ < 5.809(Q2 _ 0.1753)3[1 - 3.0098(Q2 _ 0.1753) + 4.5792(Q2 _ 0.1753)3)] 
[-13.914 + 27.679Q2 _ 13.393Q4 + 3.1924Q6 _ 0.10487^® 


if < 0.1753, 

if 0.1753 < < 0.823, 

if > 0.823, 


+ C 3 ,. • < 6.2845(Q2 _ 0.1676)3[1 - 2.9595(Q2 _ 0.1676) + 4.3355(Q2 _ 0.1676)3] 
[-15.411 + 32.088Q2 _ 17.666Q^ + 4.9355Q3 - 0.37498Q3 


if < 0.1676, 

if 0.1676 < < 0.823, 

if > 0.823, 


+ Ck^ 


2Q5 

0 


if > {Mk* +mKY, 
if < {Mk* + rriKY. 


(7) 


Let us now turn to the case of TAUOLA CLEO isospin intricate. It differs from 
the previous for the 7r“7r“7r+ only. Formulae (jS]) and (E]) are still used for but they 

are replaced by ([8]) and ([9]) in the three-prong channel. Those are also extracted from the 
TAUOLA code, but from part which was not available for the general use. 

3 Such choice is dictated by absence of scalar (like 7r(1300)) and vector (like w) intermediate states in 
this model. 

® Not all numerical values in the formula below are coming from fits. Explicitly coded numerical 
constants in TAUOLA have clear physical meaning: 0.1753=(3TO,r)2, 0.823=(Mp -|- m^r)^, 0.1676=(2m,rO + 

TO,r)2. 
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Fi = 


Ml 


Mi 


- Ml - ■ WGA{Q^) Q2 _ Ml - ■ WGAm 

Pi ■ BW{si,Mp,Tp,m 2 ,m 3 A) + P2 ■ Mp-, Tp/, m2, m3,1) 


-/33 


-Pa- 


(S3 - m|) - (si - mf) T- _ 

--- BW{s2,Mp,Tp,mz,mi,l) 


■ Pb 


(53 - ml) - (ai - ml) 

3 

(s2 - mf) - (s3 - m§) 


• BW(s2,Mp,,Tp>,m3,mi, 1) 

• BVF(si,M/2,r/2,m2,m3,2) 


■ Pb 


(Q^ + S2 — m|)( 2 m§ + 2 mf — S2) 
2 


BVF(s2,Af/2,r/2,m3,mi,2) 


- P& - -■ BW{s 2,M„,T^,1713, mi, Q) 

+ /37 • ^ • SVF(s2,M/o,r/p,m3,mi,0)y 


( 8 ) 


Again, F 2 has the same functional form as Fi. The only difference is interchange of indices 
1 and 2 for its arguments in eq. o. and that constant C 2 has opposite sign to ci. 


F 3 = 




Mi 


- Ml - • WGAm Q 2 _ Ml - ■ WGAiQl 

. (52 - mi) - (53 - mi) ^ 

_ . (53 - mi) - (51 - ml) _ p^^ 


(s2 - ml - (s 3 - m§) 


BTT(si,Mp/,rp/,m2,m3,1) 


^ (»3 '"a) ^ (»i "»?) . („, M,., r„,, m3, m,, 1 ) 


- ds 


18 s 2 

(Q^ + S2 - m2)(2m| + 2 mf - S2) 

18^^ 

-2 


BW{s 2, Mf^,Tf^,m3,mi,2) 


+ Pe- -^ ■ BW{si,Mcr,T„,m2,m3,0) 

+ de ■ ^ ■ BW{s2,Ma,Ta,m3,mi,0) 

- Pt ■ ^iT(si,M/o,r/p,m2,m3,0) 

- d? ■ ^ • BTT(s2,M/o,r/p,m3,mi,0)y 


( 9 ) 
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One can ask, why the isospin symmetry was imposed taking into account hnal states 
(that means distributions only) in the version of the code used as default in TAUOLA for so 
many years, even though it is not the property of the model used, see eg. Ref. [H]. We will 
return to this point later at a time of discussing numerical results. Current of eq. ([5]) and 
(1^ is used in TAUOLA for one-prong and three-prong of the three-pion channels in most of 
application^. Note, that neither Ref. m nor 110 ] give functional form of r —)■ 27r vr+z/,- 
current, also Ref. [ 6 ] presents functional form of hadronic current for r —)■ 2Ti^n~Uj. only. 
Parameters used in currents, will be presented later. 

5 Analytic form of hadronic current of TAUOLA BaBar 

This model became public recently, it was developed by the experimental collaboration. 
Mentioned in Sectiond] isospin symmetry and limitation of quality of the data was probably 
behind the choice of form factors used by BaBar as well. They are the same for r —)■ 
27i^7i~h'r and r —)■ 27i~7i~^i'r- BaBar model does not include resonances /o, /2 and a. Note, 
that BaBar choice was introduced after the CLEO currents (both variants) were known and 
available for them. Equations in this section are extracted from TAUOLA code numerically 
equivalent to the one used by BaBaild]. 


Fi = 


Mn 

ai 


02 _ M 2 -iM GFUNjQ^) 

'V -'*'^01 ‘-'"ai gFUN(MI^) 

{BW{si,Mp,Tp,m2,rn3,l) + j^BtP(si,Mp/,rp/,m2,m3,1) if si > (m,^ +TO,r)^, 
' 1 _I_ §1 _ W' if Si < Cm -\-m 

( M2-S1 ^ i+/3i M'p-si > 


( 10 ) 


F 2 {s 2 ) coincide with Fi(si). It has an opposite sign as well; = 0; F 4 = 0; F 5 = 0. 


GFUN{x) = 


4.1(a; - - 3.3(a; - 9mlJ] + 5.8(a: - if a: < {Mp + 

x(1.623- 


TTq / L * V TTo / 

10.38 _ 9^ I 0.65 \ 

X ^ 


TTq) j , 

if a; > (Mp -|- m.^-)^. 


Note, that eqs. CD and (I7|) aim at Oi phase-space dependence, but eq. (HD does so 
in simpler manner. CLEO model uses more sophisticated polynomial interpolation with 
inclusion of K*K production threshold; above ai can decay relatively abundantly into pair 
of kaons and a pion. 

Note, that later the BaBar collaboration was using both of these options as alternative to their default 
parameterization. 

The authors would like to acknowledge Swagato Banerjee and Tomasz Przedzinski for cross-validation 
of our modifications of the r —>■ inv hadronic current called TAUOLA BaBar, with the TAUOLA installation 
in BaBar framework at a statistical level of lOOM events. For comparisons standard distributions prepared 
with MC-TESTER [19] were used. 
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6 Resonance Chiral Lagrangian currents 

Our first hadronic current based on the Resonance Chiral Lagrangian scheme - TAUOLA 
RChL 2012 was prepared for Ref. I2D1- Similarly to parameterization of BaBar, the same 
form of current was used for one- and three-prong r decay modes into three pions. As an 
input for hts, only two, one-dimensional histograms of BaBar data were used. Later on, 
the invariant mass of n~7r~ pair became available. It turned out that contributions which, 
under isospin symmetry, provide different form of the current for one- and three-prong 
decay modes was useful. We refer to this later parameterization as TAUOLA RChL; it 
is documented in Ref. [Hj. 

We will return to isospin symmetry aspect, when we will discuss sensitivity of different 
distributions to specihc parts of current contributions. This may hint, why isoscalar terms 
introduced by CLEO many years ago, were actually dropped out in currents, such as 

TAUOLA BaBar and TAUOLA RChL 2012. 

7 Practical aspects of model construction 

There are two main paths to be taken in model building. Empirical approach is to prepare 
something describing the data in the best possible way, without considering how physical 
process works. Other choice is to start from theoretical principles, obtain distributions 
and introduce only necessary adaptations at a time of comparison with data. Of course, 
merging those two approaches is desirable and even essential. There are unquestioned rules 
like Lorentz invariance which have to be obeyed. All models used in TAUOLA and in our 
paper, respect such rules. 

In particle physics we expect Breit-Wigner distributions to describe shapes of exper¬ 
imental distributions. One should bear in mind that precision of experimental data of 
today can be very high. Multiple millions of events were collected for individual r decay 
channels. On the other hand, once inspection of the assumptions behind theoretical models 
is performed, rather modest estimations for systematic uncertainties of the predictions are 
obtained IS]. For example, for RChL models, in Ref. 1201 Section 7.4 it was argued that on 
the basis of theoretical arguments alone, the predicted uncertainty can in principle be as 
high as 30 %. Even behind such a basic principle like isospin symmetry uncertainty at the 
level of 5-10 % can be expected at least in some cases, see, e.g.. Ref. [T] Section 5.11. That 
is again more than one order of magnitude worse than precision of the experimental data. 
Comparing to differences in results of Section 19.21 it is about the same order of magnitude. 

That is why, comparison of model’s results with the data can be an exciting source 
of inspiration for further model developments, even if attained agreement with the data 
is not perfect at the beginning. Such development can be particularly promising for the 
models precisely representing distributions of several r decay channels and/or hnal states 
in e’''e“ —?• hadrons. 

In case of TAUOLA CLEO, at start, model of Ref. [22] was used. With time, as 
a consequence of comparisons with the data, it evolved into form where Breit-Wigner 
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resonances seen in experimental data were added. Some of the model assnmptions, in 
particular numerical values of the chiral couplings in soft limit were compromised. One 
should keep in mind that only at very low energies chiral dynamics is expected to determine 
unambiguously the hadron form factors probed by semileptonic r decays [231 El Eg. 

At present, more elegant approach is taken; the RChL model. It was started from 
theoretical calculations accounting to expansion in number of colors Nc ^ oo [26l \T7\ . 
Later, this model required empirical correction by the a meson introduced in Ref. [ll as 
an additional Breit-Wigner. Authors are still working out a way to include a and possibly 
other isoscalars within RChL scheme using unitarity and analiticity constraints to include 
this re-scattering effect. 

Other interesting approach in modeling was taken, e.g., in Ref. [28]. It describes 
hadronic current for r —)■ This modeling is in part data driven. As measurements of 

e’''e“ annihilation into four pions arrived [29], they could be used to model distribution of 
invariant mass of the whole dvr system^. Authors introduced correction into matrix ele¬ 
ment to reproduce experimental distribution. Then, aivr dominance as intermediate state 
was assumed. From that point on, any desired modeling of ai could be used. The one of 
Ref. [30] was chosen as it was applied for data as well. 

Those are only a few of many possible approaches. There is no straightforward way of 
telling that one is better than other. Actual choice has to be motivated by quality and 
precision of experimental data. Let us address this point in the next section. 

8 Data representations and numerical methods 

Experimentalists and theoreticians may prefer distinct forms for data representation in 
model building, comparisons and hts. In this section we will concentrate on methods we 
use and those that were used for constraining investigated models in the past. 

The most simple way to compare models and also to ht the experimental data is to look 
only at the total width of the channelil (we will call it option 1). Such approach is used 
when there is no sufficient data available for a more advanced fit. Even though it is only 
one number, it can be used as a first test of particular model. Total width is usually well 
documented by PDG tables, see, e.g.. Ref. im including careful analysis of experimental 
statistical and systematical error. 

Next approach is to look at one-dimensional distributions. In this case, one should seek 
the possibility of having experimental data for all possible invariant mass distributions. 
This is not always possible. Using as an example the channel, one must be aware 

that it is easy to have many models that fit perfectly, let’s say the 7r“7r“7r+ spectra, while 
having distributions of and off. We can articulate 3 options of one-dimensional 

distributions in use for the case of r —)• dvru,- current input: 

It was done by relating 27r+27r“ and 7r+7r“27r° systems produced in e+e“ annihilation with 7r+27r“7r° 
and 7r“37r° produced in r decays. For that, isospin symmetry was applied to distributions measured from 
e“'"e“ annihilation giving prediction on shapes of four-pion system invariant mass spectra in r decays. 

Option 1 was used, e.g., in Ref. m for Btt currents. 
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• system invariant mass onl}0 (option 2). 

• 7r“7r“7r+ and 'k~'k^ invariant massed (option 3). 

• 7r“7r“7r’'‘, 'k~'k^ and n~TT~ invariant massed (option 4). 

In more advanced approach of Dalitz distributions (option 5) slices in were 
used already in Ref. [6]. This way of data representation is in nature three-dimensional^. 
Models originating from the CLEO collaboration were htted in that way. To explore the 
whole structure of the decay channel methods of Ref. [4] have to be introduced (option 
6). Up to now, no complete analysis was performed with their help. Only Ref. [6] mentions 
cross-check with such methods. 

Let us turn to practical aspects of our paper comparisons for distributions with options 
listed above. For that purpose we describe re-weighting method. We use it often due to its 
numerical efficiency and simplicity. When generating events from Monte Carlo generator 
for different models, statistically independent event samples are produced. But one may ask 
what is the probability of getting particular event with certain set of momenta depending 
on the model we use. Having all momenta of particles in an event, one can calculate 
the amplitude for such specihc conhguration. Therefore, one can calculate amplitudes for 
the event using each model and compare. This leads to the possibility of re-weighting 
events by using one model as reference and instead of generating a new sample with a 
different model, calculating amplitudes from that other model for existing events. Ratio 
of new amplitude squared to old one, can be then attributed to the event as its weight. 
Thanks to weights use, samples are correlated. Statistical error affects the difference for 
the compared models only. Using this strategy for a different model we may not only get 
the new distributions with different shape but also its integral, that is the total width of 
the channel. Method can be applied to minuscule changes of model parameters, which are 
needed, e.g., for calculation of derivatives. 

Re-weighting events can be used to estimate the contribution from parts of the model 
amplitude as well. The ratio of investigated part of the amplitude squared to a whole 
amplitude squared provides then the weight. Evaluation of interferences between parts of 
currents is also possible on the event by event basis. One has to evaluate how far from 1 
is the sum of weights coming from all parts. 

Let us look into that in more detail. For any integrated distributions bin always repre¬ 
sents average weight of events hlling it. Once histogram is hlled, information from separate 
events is lost. Because of that, interference observed in any histogram is reduced with re¬ 
spect to its full size. Effect in integrated distribution is dependent on selection of bins. 
By comparing interference coming from the event by event analysis and integrated distri- 

Similar to option 2 fit was performed for 47r current in Ref. |28] with only invariant mass of four 
pions used. 

Option 3 was used in fits of TAUOLA RChL 2012. Probably also BaBar model was fitted in this 

way. 

Option 4 was used for fits of TAUOLA RChL. 

There are 8 plots, each in different range of on x-axis there is Si and S 2 is on y-axis. si is always 
chosen as the bigger of the two possible 7r“7r+ pair invariant masses. See plots of appendix 0 
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butions one may evaluate how viable used distributions are to constrain interference and 
if model can be explored with available data sufficiently well. Integrated contributions 
from parts of current are identical for any distribution, but interferences may differ vastly. 
Analysis of option 1 is expected to give very poor insight on interferences, and options 2-4 
can give only limited information. Option 5 is expected to give the best way to explore 
models in that regard. Yet, it still does not explore the whole structure of multidimensional 
phase-space, contrary to methods of Ref. [1]. 

9 TAUOLA CLEO, TAUOLA CLEO isospin intricate, 
CLEO publ. 

Beside functional form of currents, numerical parameters affect shapes of distributions. 
Among them, mainly masses, widths, and couplings of the intermediate resonances differ. 
That is why, we collect all constants to the last detail, taking into account variants in the 
implementations. The set of parameter^ is extracted from the TAUOLA code, to make 
sure that they were used in TAUOLA CLEO and TAUOLA CLEO isospin intricate 
currents. These currents rely on analysis of data using option 5, but further cross-checks 
using option 6, were also mentioned in CLEO publications. Published experimental data 
for T —)• were used by the collaboration. Predictions for r —)• were 

documented only in conference contributions [TOl dU. It is not the theoreticians or MC 
authors right to decide whether predictions for the 7r“7r“7r''' channel should be obtained 
from the case at the level of distribution or taking into account isospin rotation 

of intermediate states. Actual choice must rely on analysis of systematic errors. Our 
default current TAUOLA CLEO implies that one considers contributions from interme¬ 
diate states: fo, f 2 , o' as effective ways to improve distributions. We use then distribution 
level isospin symmetry to get predictions for 7r“7r“7r+ from we use the same cur¬ 

rent for the two cases. Second variant, TAUOLA isospin intricate, was also prepared 
by the collaboration but was not published. We have nonetheless distributed it together 
with all versions of TAUOLA, but not made available for direct TAUOLA usJ^. To activate 
it, modification of code was always necessary. Let us stress that the numerical difference 
between the two options is comparable/smaller than the differences to other parameteriza- 
tions presented in later sections. We believe also that the choice lies within the systematic 
uncertainties. This would of course have changed with CLEO publication or with the fu¬ 
ture publications of other collaborations. For further discussion and numerical comparison 
of different current options see the following subsections and Sections dH [T^l 

i®Note: we present here only those parameters that are different as only those can introduce any dis¬ 
crepancies between CLEO publ. and TAUOLA CLEO. The complete list of parameters is given in 
Section [TTl 

There are other options of hadronic currents distributed, but never made available for general use 
as well. The aim was to stimulate future work or for archivization. That was for example the case of 
Tr~Tr~TT'^TT^ current used for background modeling in Delphi analysis of Ref. [32] and (Atttt)" currents of 
OPAL [33]. Alternative ALEPH parameterization of hadronic currents was documented in Refs. 019]. 
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The model described in Ref. [ 6 ] and the implementation used in TAUOLA differ only by 
parameters stored in Table [T] One should stress that Ref. [ 6 ] is for the mode, 

whereas TAUOLA CLEO use the same current for the mode as well. 


TAUOLA CLEO 

GLEO of Ref. [ 6 ] 

Mp = MRO =0.7743 GeV 
Fp = GRO =0.1491 GeV 
Mr = 1.777 GeV 

Ma, = 1.275 GeV 

F^^ = 0.700 GeV 

Mp = MRO =0.774 GeV 

Fp = GRO =0.149 GeV 
not mentioned in the paper 
Ma, = 1.331 GeV 

F„^ = 0.814 GeV 


Table 1: Parameters differing in TAUOLA and Ref. [6]. These differences in parameters are 
the reason of differences visible in Figure [H 

Analytic form of TAUOLA CLEO isospin intricate current for r —)■ can 

be understood as a simple consequence of isospin symmetry applied to the r —)■ 
case. But we must make sure if the hadronic current to be transformed by isospin rotation 
is sufficiently established, see discussion in section 114.11 

One could wonder, if it is not better to drop the terms of a, /o and /2 from the current 
parameterization, so currents for both channels would have the same functional form. 
Indeed, that was often the case. For example, parametrizations TAUOLA BaBar and 
TAUOLA RChL 2012 introduced later, used identical currents for and . 

In TAUOLA CLEO the same current for both channels was used disregarding the fact 
that it was seemingly introducing a of a charge ±2. In fact, it was simply assumed that 
a was just a feature of the current parametrization and not a real physical state. This is 
supported by the observation that on Dalitz plot from Fig. 3 of Ref. ra there seems to 
be a resonance-like structure in S 3 , the invariant mass of the pair. In corresponding 

plot. Fig. 5 of Ref. [ 6 ] for case, the same structure can be seen in It is most 

likely originating in part from interferences, but it can be parameterized as a resonanc^. 
Let us now look into numerical differences between basic variants. 

9.1 Numerical differences of TAUOLA CLEO and CLEO publ. 

The MC sample presented in [B] is not available now, therefore we can only try to recreate 
results. To do so, we need to identify all differences in initializations. It appears, that only 
a few parameters differ, see Tabled] Resulting difference is small but clear on, see Fig. dJ 
Note, that these two variants used hxed values for masses and widths of resonances (while 
RChL and BaBar models htted them). We can only speculate it is either due to stabilitj0 
of fit procedure or simply the belief that use of results from other (non r) measurements 
is better. Apart from extreme kinematic regions, the largest differences can be observed 
in 7 r“ 7 r“ 7 r''' distribution. 

See also Section 114.21 where we discuss some possible origins of unexpected structures. 

Overpopulation of parameters may affect stability. If that was not the case, we may suspect lack of 
computing power, as additional parameters increase CPU time needed for fit to converge substantially. 
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Figure 1: Invariant mass, GeV units, of (from the left) vr^vr’'', in TAUOLA 

CLEO with numerical constants shifted to values of Ref. [6] (red) and TAUOLA CLEO 
with default parameterization (green). For black line which represent ratio of red plot to 
green one, left hand side scale should be used. Right hand side scale represent number of 
events per bins of 10 MeV. 


We can also investigate differences in terms of three-dimensional plots. Fig, [represents 
Dalitz plots obtained from TAUOLA CLEO. As we do not have data nor MC generated 
sample from the CLEO collaboration, only human eye comparison is possible with their 
published results. Mentioned Fig. [14] of our paper should coincidS with Fig. 5 of. Ref. |6]. 
No convincing differences can be seen. 


9.2 Numerical differences of TAUOLA CLEO and TAUOLA CLEO 
isospin intricate 

These two currents have the same values of parameters but differ by functional form of 
hadronic current for r —)■ TAUOLA CLEO is represented by equations (jSjl. 

(j6]), while TAUOLA CLEO isospin intricate by eqs. (E]), (E]). The differences are large, 
see Fig. E] judging with today standards for measurements, but not larger than of relatively 
recent parameterizations TAUOLA RChL and TAUOLA RChl 2012 shown later in 
Fig. [6] In fact they are at the 10% level, that is at the level of the theoretical uncertainty 
for use of isospin symmetry [1]. 




Figure 2: Invariant mass of (from the left) ^ in TAUOLA CLEO 

(red) and TAUOLA CLEO isospin intricate (green). For details of histograms dehni- 
tion see caption of Fig. [TJ 

Note, results of Ref. |6] measurements are for the mode, but the same functional form of 

hadronic current is used for the 7r“7r“7r+ in TAUOLA CLEO, that is why such comparison is valid. 


14 














































10 Further variants of currents used by CLEO 

Currents described in this paper are not the only ones developed by experimental col¬ 
laborations. We would like to hint why certain solutions were used and documented in 
publications while others seem to be forgotten. 

We should start the discussion with CLEO publ. This current is fully described in 
Ref. [6]. But the very same publication mentions a couple of other different variants, 
stating for some of them even better agreement with the data. So why is that any of those 
did not become main model of the collaboration? 

First, we should mention variant of model where a mass and width were £t parameters. 
Best agreement^ with the data was obtained with M„ = 555 MeV and Fg- = 540 MeV 
(page 22 of Ref. |6]). Such fit caused shift in all couplings not exceeding 20%. We may 
again only speculate that it was believed that it is better to use theoretically predicted or 
experimentally well established values. Fitting of sigma mass and width was probably per¬ 
formed in order to check how seriously this resonance should be treated. Result: different 
values than used for nominal model (Mg. = 860 MeV and Fg. = 880 MeV), was most likely 
not very encouraging. 

Secondly, there is mentioned in Ref. |6] a study on pseudo-scalar contribution from 
7r(1300). It also improved the £t result, but with available data its existence could not be 
established with unquestioned signihcance. The very same reasoning was present in study 
of a'l contribution, which also was not included in nominal model. 

PhD thesis [12] presents similar study of pseudo-scalar 7r(1300) in r —)■ 
decay. Different ways of inclusion of scalar current were tested there. Similarly agreement 
with the data got better, but was not considered signihcant enough. In this study hadronic 
current used for non-scalar part was the same as TAUOLA CLEO isospin intricate. 
It suggests again, that at some point the current without pseudo-scalar contribution was 
used by the collaboration. 

Another thing under consideration by CLEO was the resonance radius. Their study 
showed small importance. Therefore, it was set to 0 in order to simplify the model. Keeping 
things as simple as possible seem to be important for many choices of collaborations. 

The BaBar collaboration had all CLEO currents available, yet decided to use simpler 
currents. Apart from simplicity, one can suspect another reason behind such a choice. 
The BaBar collaboration use commonly one-dimensional histograms in analysis even in 
our times [3l]- Inclusion of any resonances other than p and p' is not to be easily exploited 
without three-dimensional distributions. On Fig. [3] one can see, that differences are not 
very compelling. They are mostly below 10%. Comparison of data and theoretical model 
on scatter-grams like Fig. [20] can show more clearly importance of other resonances, but it 
was never performed. On the last hgure, ratio of Dalitz plots coming from BaBar model 
and CLEO model is shown. While ratio of invariant mass distributions was mostly between 
0.9-1.1, ratio in some areas of Dalitz plots escapes 0.5-2 range. 

^^Negative doubled logarithmic likelihood was 43 units below that of nominal fit. 
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11 Comparison of TAUOLA CLEO, TAUOLA CLEO 
isospin intricate and TAUOLA BaBar 

In this section comparison of CLEO and BaBar model used in TAUOLA is shown. Parameters 
are divided into three groups depending on the way they are introduced into the code and 
how they were used. Tab. [2] presents those dehned globally, outside form factors. Table 
m presents those introduced inside hadronic current only. Parameters in this table were 
kept constant during the fitting procedure of CLEO, but were fitted for BaBar modeling. 
Table m stores parameters fitted to the data. They are also defined inside hadronic current. 
Note, that some parameters can be different for diverse parts of code, for example mass of 
p resonance. This may be misleadingly worrisome, but it should not be considered as an 
error. Parameters used outside hadronic form factors affect optimization of efficiency of 
phase-space generators, that is essentially speed of generation only, or are well established 
constants like Their variation does not affect output distributions, or variations are 
simply unphysical. 


TAUOLA CLEO 

TAUOLA BaBar 

13 = BET =0 

ttIttO = AMPIZ =0.134996 GeV 
= AMPI =0.139570 GeV 

Mp = AMRO =0.775 GeV (used in GFUN) 
Ma, = AMAl =1.251 GeV 

P„, = GAMAl =0.599 GeV 
/3 = BET =0 


Table 2: Parameters dehned outside form factors (stored in INIMAS). 


The TAUOLA CLEO code is essentially as described in Ref. [6] and named ’nominal 
ht’. It differs by some numerical constants only, see Tabled! The CLEO collaboration did 
not ht masses and widths of resonances but only coupling constants (3i. Such choice was 
driven by theoretical assumptions and measurements of masses elsewhere. This may be 
the best solution when we lack good enough measurements of state used in model. Let us 
recall. Section dQi how insensitive the model can be to the a meson properties. 

11.1 Numerical comparison of TAUOLA CLEO and TAUOLA 
BaBar 

Here, we present numerical comparison of TAUOLA CLEO and TAUOLA BaBar. On 
Fig. [3l distributions of invariant masses of are shown. Looking at 

the plots one notice that except the lack of some resonances mentioned earlier, also lack 
of K*K threshold in BaBar model is visible in distribution of mass (and even 

more visible on ratio plot) in 1.4 GeV region. Moreover, diherent Oi mass plays a role. 
In addition, on Fig. dOl of Appendix O ratio of Dalitz plots is shown. Lack of /o, / 2 , cr 
resonances in model used by BaBar can be also seen there. Diherences are higher than for 
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TAUOLA CLEO 

TAUOLA BaBar 

Mp = MRO =0.7743 GeV 

Mp = ROM =0.773 GeV (used in Breit-Wigner) 

Tp = GRO =0.1491 GeV 

Pp = ROG =0.145 GeV 

Mp, = MRP =1.370 GeV 

Mp, = ROMl =1.370 GeV 

Tp, = GRP =0.386 GeV 

Tp, = ROGl =0.510 GeV 

Mf^ = MF2 =1.275 GeV 

Tf^ = GF2 =0.185 GeV 

Mf^ = MFO =1.186 GeV 

P/o = GFO =0.350 GeV 

M„ = MSG =0.860 GeV 
p^ = GSG =0.880 GeV 
rriT,^ = MPIZ =0.134976 GeV 
= MPIG =0.139570 GeV 
rrixo = MK =0.496 GeV 
= MKST =0.894 GeV 
Ma, = AMAl =1.275 GeV 
Pai = GAMAl =0.700 GeV 
= BTI = 1 

(3i = BETAl =-0.145 


Table 3: Parameters used in definition of hadronic currents in TAUOLA; both variants of 
CLEO parameterization and BaBar parameterization. This table stores parameters that 
were not fitted by the CLEO collaboration, but were taken from PDG tables or theoretical 
predictions. Corresponding BaBar parameters are stored here for completeness. 


comparison of TAUOLA CLEO isospin intricate and TAUOLA CLEO, see Fig. [T] 

Invariant mass of system exhibits biggest difference in terms of ratio. 





Figure 3: Invariant mass of (from the left) tt tt’*', tt tt , tt tt tt’*' in TAUOLA CLEO 
(red) and TAUOLA BaBar (green). For details of histograms definition see caption of 

Fig.E 
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/32 = BT2 = 

/^g = BT3 = 
j3^ = BT4 = 0.87e*°-53^ 

(3^ = BT5 = 0.71e*°-5®^ 

/3g = BT6 = 2.10e*°-23^ 

(3^ = bT 7 = 0.77e-*°-5^’^ 
C-g^ = C3PI = 0.23842 
Ck* = CKST = 4.76212^3^ 


Table 4: Parameters used in definition of CLEO hadronic current in TAUOLA obtained by 
the CLEO collaboration from fit to the data. 

11.2 Numerical comparison of TAUOLA CLEO isospin intricate 
and TAUOLA BaBar 

For completeness we provide comparison of TAUOLA CLEO isospin intricate and 
TAUOLA BaBar. Results are very similar to previous comparison and can be seen on 
Fig. m We do not discuss those results as no new aspects appear. Largest difference is 
present in invariant mass spectrum. 



Figure 4: Invariant mass of (from the left) , 7t~7t~, 7r~n~7i~^ in TAUOLA CLEO 
isospin intricate (red) and TAUOLA BaBar (green). For details of histograms defini¬ 
tion see caption of Fig. [TJ 


12 Comparison of TAUOLA CLEO and TAUOLA RChL 

In this section, numerical comparison of CLEO and RChL model used in TAUOLA is shown. 

We will not recall analytic form of RChL current. Full description of TAUOLA RChL 
2012 current is given in Ref. [20] (see also [I3]), and for TAUOLA RChL in Ref. [14] . 
CLEO and RChL models differ a lot in a way of construction and also in parameters, but 
one may easily guess what affects the spectra in a way as seen on Fig. [S] Most importantly, 
inclusion and parameters of a resonance are different. The a of lower mass and width 
than in the CLEO model, affects mainly low mass region of in RChL, while 

in CLEO model it affects the whole spectra. Second, clearly visible difference appears 
in 7r“7r“7r’'' invariant mass distribution. It comes from different modeling of ai. CLEO 
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translated theoretical modeling into polynomial of eq. ([7]). RChL keeps direct model 
coded, bnt parameters of K* meson were not fitted, as the 7r“7r“7r’'' channel shonld not 
be very sensitive to them. Reasonable modihcation of those parameters wonld reqnire 
simnltaneons £t of the KKn channel. Snch fit conld possibly rednce difference in area of 
the K*K threshold mentioned in previons section. 


( A 

-,5*1 
-i 220 

-;200 3 

-il80 

-il60 2.S 

-il40 


-^70000 2.4 

] 2.2 

H60000 2 

H50000 

4 1.6 

H40000 1-4 

A 1 

1 \ : 


-i120 

-ilOO 1.5 

-jso 


3 1.2 

H30000 1 

3 O.B 

-320000 



-!«> 0.5 

-i20 

-q 10000 

] 0.2 

do 0 

i..J. . Aj 


Fignre 5: Invariant mass of (from the left) tt tt’*', tt tt , tt tt tt’*' in TAUOLA CLEO 
(red) and TAUOLA RChL (green). For details of histograms definition see caption of 

Fig. [H 

In Fig. HI] of appendix [U] we can hnd three-dimensional ratio of RChL to CLEO model 
resnlts. On that plot, we can see areas where difference between the models exceeds 
factor of 2. This is mnch bigger difference than seen on one-dimensional distribntions. As 
CLEO model was htted to three-dimensional distribntion it is a very worrisome observation 
from the RChL perspective. It may hint that even thongh one-dimensional BaBar data 
distribntions are described well na, nnderlying physics model may be imperfect. Actnal 
jndgment shonld not be completed withont comparison with experimental distribntions of 
higher dimensionality. 


12.1 Comparison of TAUOLA RChL 2012 and TAUOLA RChL 

For completeness, we present comparison of hrst version of RChL model and a new one. Be- 
canse of the comparison with the data, especially inability to obtain satisfactory agreement 
in invariant mass, when all invariant mass distribntions were htted simnltaneonsly 

TAUOLA RChL 2012 was modihed. Those two models differ mainly by inclnsion of a 
resonance in the TAUOLA RChL version. Other differences are model parameters and 
minor correction for Conlomb interaction. Resnlting differences in distribntions can be 
seen on Fig. [6l 

At this point we mnst emphasize that those two versions of RChL model were htted 
to diherent data-sets. Invariant mass of system became available for hts only after 

development of TAUOLA RChL 2012. This model was htted to invariant masses of 
Tr^TT'*', 7r“7r“7r’'‘ systems only, while TAUOLA RChL to invariant masses of , 7r~7i~, 
7T~7r~7T~^. It was not the only diherence. Also binning of histograms changed. New version 
was htted to histograms of lOMeV bin widtf@. Distribntions nsed in ht of TAUOLA 
RChL 2012 had bin width of 20MeV. 

This kind of binning is used in all one-dimensional plots of our paper. 
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Figure 6: Invariant mass of (from the left) 7r“7r+, 7r~n~, in TAUOLA RChL 

(red) and TAUOLA RChL 2012 (green). For details of histograms dehnition see caption 

of Fig. [H 

13 Numerical effects of interferences within currents 


So far, we have concentrated on presentation of different currents used for Monte Carlo 
simulations of r decays. We have looked into their analytical form as well as the numerical 
values of the parameters fixed by confrontations with experimental data. Comparisons 
obtained from different parameterizations were reported as well. In the present Section let 
us look into numerical consequences of some features of the currents. We will concentrate 
on interferences between distinct parts of the current, which can be separated, with the 
help of isospin transformation. 

Final states of long-living scalars with given multiplicity and charge combination, can be 
formed through intermediate resonances used in description of hadronic current. One can 
not think of them as sum of non coherent contributions. Interferences represent important 
aspect of hadronic current which should be understood and controlled. Even if used for hts 
experimental distributions are correctly reproduced, unnatural response to modifications 
of individual contributions to the current and interference can hint for modeling flaws. 

In principle, addition of a new resonance to the hadronic current should result in neg¬ 
ative interference below the peak and positive abovj^. There is no clear prediction on 
how strong such resonance should interfere. In case of a narrow intermediate resonance we 
expect such interference to cancel out in contributions to total rates or relatively inclusive 
distributions. Especially in case of wide resonances, we may model a process with some¬ 
thing that is not of proper dynamical nature but still describes data quite well, simply by 
accident. 

In this section, we will investigate RChL and CLEO modeling of the three-pion channels. 
We concentrate on interference of part depending on isoscalar^ with the remaining part 
of the current. These two parts of hadronic current transform differently under isospin 
rotation. In principle it is also a test on how justihed (or not) is isospin rotation from 
7r°7r°7r“ to currents. 

Let us look at the contribution to the rate. For the p part, the cr, / 2 , /o part 

That is true in case when resonances come in with the same sign. In our case signs are opposite. 

Contribution only from a in RChL, while for CLEO from u, / 2 , /o- 


20 
























and their interference contribntion are respectivel}@ 75.1%, 4.4% and 20.5%. Usage of 
TAUOLA CLEO isospin intricate cnrrent for 7r“7r“7r+ gives 74%, 5.8% and 20.2%. 
Interference looks slightly different for fnlly differential distribntions. Then, it is 22.6% 
in TAUOLA CLEO, and 20.7% for TAUOLA CLEO isospin intricate case, when 
method as explained in Section [H] is nsed. For the three-dimensional scatter-grams of the 
CLEO collaboration style, see AppendixOfor details, the interference is respectively 22.3% 
and 20.6%. This demonstrates nearly perfect sensitivity of snch a method. 







Fignre 7: Invariant mass of (from the left) 7i~7r~^, 7r“7r“, nsing TAUOLA CLEO 

isospin intricate (top) and TAUOLA CLEO (bottom). Green plot represents contribn¬ 
tion from isoscalars {a, / 2 , /o) and bine from remaining part. Red plot represents the whole 
cnrrent. Bine and green plots wonld snm np to red one, in the absence of interference. 
Percentages of the three contribntions are given in Table O 


Current 

Absolute width 

ID distributions 

3D distributions 

Absolute width differential 

ir^Ti-^TT- CLEO 

75.1% 4.4% 20.5% 

75.1% 4.4% 20.8% 

75.1% 4.4% 22.4% 

75.1% 4.4% 22.7% 

7r~7r~7r~^ CLEO 

75.3% 4.3% 20.4% 

75.3% 4.3% 20.7% 

75.3% 4.3% 22.3% 

75.3% 4.3% 22.6% 

7r~7r~7r’*” CLEO isospin intricate 

74% 5.8% 20.2% 

74% 5.8% 20.2% 

74% 5.8% 20.6% 

74% 5.8% 20.7% 

7r~7r~7r'^ R,ChL 

97.5% 6.5% 4% 

97.5% 6.5% 11% 

97.5% 6.5% 18.2% 

97.5% 6.5% 18.5% 

ir“7r“7r- RChL 

94% 1.3% 4.7% 

94% 1.3% 7.4% 

94% 1.3% 9.5% 

94% 1.3% 10% 


Table 5: Comparison of contribntions from the isoscalars from the remaining part of the 
cnrrent and from interference, depending on a way of investigation. In each cell, hrst 
nnmber is from the spin one intermediate states, second is from isoscalars and the third one 
is interference. For 1-D (one-dimensional) distribntions interference effect of the histogram 
with largest effect is taken. In the case of ” Absolnte width differential”, average of (module 
of) interference over all events is shown. 

In Table [5] we collect such results also for other versions of the currents. One must 
remember that the interference visible in histograms is somewhat dependent on binning. 
We do not have the same binning as the CLEO collaboration; they had less bins per energy 

When the same current is used for 7r“7r“7r+ corresponding contributions are 75.3%, 4.3% and 20.4%. 
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unit in every distributions. Therefore, while they probably could have seen the same 
behavior, the interference observed was most likely smaller. Interference is an important 
numerical property of the distributions which can not be ignored while making the choices 
for phenomenological approach^. 

Having discussed CLEO modeling let us investigate the same problem in RChL model. 
On Fig. [HI one can compare contributions from different parts of current on one-dimensional 
distributions. We can see that plot follows theoretical expectation on how interference 
should look. Also, there is much less interference than for CLEO models, see Table [5l 
Moreover, size of interference for the one-prong and three-prong channels is not as similar 
as in CLEO modeling. While both CLEO models have almost the same level of interfer¬ 
ence in both channels, RChL has twice more in the three-prong channel. Also, size of a 
contribution differ a lot between the one- and the three-prong channels, while it is close to 
the same in both CLEO modelings. Another property worth noting is limited sensitivity 
of RChL model to interference when only one-dimensional distributions are used. It means 
that the most recent ht of this model does not fully control^ influence of a and fit to 
three-dimensional distributions is strongly desired. 

As both models aim at describing the same decay channels we should be concerned 
that interplay of internal parts of the current is quite different. It proves that we do not 
have good understanding of physics behind it, especially what exactly is the a state, in at 
least one of the two approaches. Both models hnd it important 0 la to include a, but 
its description is far from consistently. 


14 Interplay of experimental and theoretical input 

In the process of creating a model for hadronic current to be used in MC there are two 
aspects. Theoreticians provide model for hadronic currents and experimentalists use it and 
fit parameters. Usually, theoreticians concentrate on consistency of formulation and prefer 
to avoid ad hoc phenomenological corrections. Very often, model building is disconnected 
from data analysis. On the other hand, experimentalists would like to have model perfectly 

Later, in particular in Section 114.21 we discuss how resonances are deformed and even can produce 
artificial peaks. Therefore, resonances are not always needed to emulate peaks present in experimental 
distribution. Also, because shapes of resulting modifications to distributions are distinct with respect to, 
e.g., Gaussian peaks of resonance parameterizations, we can be concerned about choice between TAUOLA 
CLEO and TAUOLA CLEO isospin intricate. This especially, as dominant contribution of isoscalars 
is a result of interference. Confrontation with the experimental data should be used to resolve this matter. 

The reason why interference effect is better visible in one-dimensional distributions in case of CLEO 
model than in RChL is technically quite simple. In the first case the mass of the a is larger so the region 
of phase-space above the peak is restricted, no cancellation occurs. Nonetheless it is important to realize 
such features while evaluating suitability of experimental input for the model fit. It should not be assumed 
that the conclusion will remain the same for all range of the fit parameters. 

We may speculate that in part, this difficulty in control contribution from isoscalars was a a reason 
behind introduced later much simpler current in BaBar Monte Carlo simulations. Note, that in Ref. m 
parameters of a contributions were featuring errors from the fit larger than for other resonances. We will 
return to this point in Section 1141 
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Figure 8: Invariant mass of (from the left) ^ 7r“7r“, 7r“7r“7r''' using TAUOLA RChL 
current. Green plot represents contribution from isoscalars and blue one from remaining 
part. Red plot represents whole current. Blue and green plots should sum up to red one in 
the absence of interference. One can see negative interference for RChL model whereas it 
was almost everywhere positive for TAUOLA CLEO isospin intricate, see Fig. [71 and 
Table [5l 

describing the data. At the same time, they must have control over any possible sources 
of experimental errors. Therefore, they may prefer to sacrifice some model assumptions 
for the sake of simplicity and more straightforward systematics analysi^. As experiment 
controls data sample it is experimentalists to give ’’stamp of approval” for the hadronic 
current. Such thing was done for 7r°7r°7r“ current of CLEO modeling, but is lacking for 
7r“7r“7r''' case. Only conference contributions describe some modelings of this channel. In 
case of the CLEO modeling we may only speculate if there were some unresolved problems 
with this channel or it was simply lack of the manpower that stopped the collaboration 
publication from appearing. In the following subsections we will discuss some aspects of 
model-building that can play a role. 

14.1 Ambiguities of experimental inputs 

Isospin symmetry between r decay channels, can be restricted to a relation of final states 
only. As a consequence, for intermediate states, such a choice could be understood as 
including contribution of e.g., a —)■ TT~n~ for the channel r —)■ . This of course 

would mean an error. However, as a was not well established (and very broad) resonance 
at the time of the CLEO work, this may not be unreasonable. Back then, /o(400-1200) [35] 
with width ranging from 400-1000 MeV was expected. Later on, in years 2002-2012 /o(600) 
was expected [36]. Recently, in 2012 it was shifted to /o(500) [37], thanks to numerous 
calculations, e.g., [38l [39] claiming both mass and width close to 500 MeV. In current 
measurements this resonance remain extremely wide, while the CLEO collaboration used 
even wider (880 MeV) and more massive (860 MeV) cr, as of Ref. [40] . 

That is why, from the experimentalist point of view, one can think of this amplitude 
contribution as a heuristic modihcation for the current, of not much quantum number 
meanin^^. Such interpretation may hold, unless it is excluded by comparisons with exper- 

In particular that there are no strongly correlated parameters in the fits. 

In Ref. [5] we read: The form factors as defined in Eqn. A2 (Equation A2 of Ref. corresponds 
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imental data. In the meantime, such an option in the Monte Carlo, can be helpful to study 
if the contribution is indeed the a. Alternatively, it is possible to think that broad a is 
just a pretext to correct current of oi, p intermediate states, thus may be identical for one- 
and three-prong r decay modes. Similarly one can think of /q. The /2 is not as broad, but 
even in this case contributions to si, S 2 spectra are not clearly localized in the phase-space. 
We should also mention that both /2 and /o have limited phase-space for decay, therefore 
their Breit-Wigner distributions are strongly deformed. On top of that, in Section [13] we 
have shown that isoscalars in the models we have presented, do not contribute as separate 
resonances to the decay. They contribute positive interference over almost whole spectra. 
That is why, such imposed isospin symmetry understood as use of the same current for the 
7r“7r“7r''' and the decay channels of r can be justihed. Let us recall that isospin 

symmetry may be broken and of limited predictive power starting from 5-10% precision 
level (Ref. [1] Section 5.11). Also, there is very small contribution to the currents from 
isospin violating decay^. 

14.2 Further concerns for model-data confrontation 

Let us look into details of the TAUOLA CLEO isospin intricate current, and the 
decay channel r —>■ 7r“7r“7r+i/^ . The contribution from isoscalars may be missed from the 
distribution of invariant mass spectrum: no clean features are seen. Fig. [9]green line. 

If we artihcially increase by a factor of 10 isoscalars amplitude, they dominate the shape 
of TT~'K^ invariant mass spectrum, see red line, but clear peaks are still absent. Position 
of the peaks is shifted with respect to the masses used in the propagatoi0 of hadronic 
current. One can see, blue line, that two peaks emerge clearly if widths of isoscalars are 
reduced by a factor of 10 as well. Looking at all three lines, one realizes that appearing 
low energy peak/deformat ion present in TAUOLA CLEO isospin intricate (in other 
modelings as well) is in a big part indirect. There are two possible combinations of 
and both are put into the histogram. Lower peak-like bump comes from small phase-space 
available for ti~ when system of other vr” and 7r+ is produced with high invariant mass. 
That forces low invariant mass bump for the 7r“7r+ system, even though those two pions 
do not originate from a common intermediate resonance. 

Resonances at the end of the phase-space represent challenging part of the description. 
As we could see from Fig. (Q] their shapes are substantially deformed from the ones of 
Breit-Wigners. The other effect is also possible. In Ref. |13] it was demonstrated how 

to our eq. O/ do not have simple correspondence with those that can be assoeiated with speeific resonant 
contributions to the hadronic current, which seem to support our statement. 

In the T —>■ decay mode isospin violating signals are expected mi- There are two decays that 

can feed Stt spectrum. For the three-prong channel it is r —>■ un:~Vr with subsequent w —>■ 7r“7r+ decay. 
The one-prong channel can be fed by the decay r —>■ r]TT^TT~Vr (due to p — mixing). Those decays are 
estimated to give small contribution to the total rate at the level of 0.4% and 10“^%, respectively. Studies 
of T —)> riTT^TT~i'T case were also done within RChL scheme [42]. 

^"^See formulae (Hj), (Hj), and Tables jUE][D for the definition of p, a, A, and /o propagators as used by 
the CLEO collaboration. 
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Figure 9: Invariant mass of system in r —)■ 7r“7r“7r+ decay using TAUOLA CLEO 
isospin intricate (green), current with multiplied 10 times amplitude from isoscalars (red) 
and also with width of isoscalars reduced 10 times (blue). 


deformations due to resonances acting as phase-space constraints, may result even in arti- 
hcial resonance-like peaks. This observation is another reason why one may be hesitant to 
exploit isospin symmetries of currents in an uncritical way. 

Another concern is the size of contribution due to interference. It may be much larger 
than the effect of extra resonances alone. At the same time, it is not fully explored by 
inclusive distributions. See Table [5l for RChL, interference effect is much larger for the 
three-dimensional histograms than for the one-dimensional ones. For the CLEO case, the 
difference was rather minor. One has to remember of such effects, while introducing new 
resonances into the currents. It must be carefully studied if they affect shapes of the 
distributions in a controlled way. 

Finally, if numerical effects are small and weakly depend on model parameters, Monte 
Carlo methods may not be best suited for the studies. Especially if systematic errors are 
to be taken into account [l^. If available, semi-analytical methods are better. Therefore, 
we have prepared special arrangement which is briefly presented in Appendix 


15 Summary 

In this paper, we have investigated several variants of parameterizations of hadronic cur¬ 
rents used in Monte Carlo programs for simulations of the r lepton decays. As a reference 
for comparisons we have used TAUOLA CLEO parameterization, because of its well 
established status; it is supported by the collaboration publication. The hrst investigated 
variant was TAUOLA CLEO isospin intricate, as prepared for TAUOLA but not com¬ 
monly used. This parameterization is documented by conference contribution only. We 
have evaluated the numerical size of the differences, which were found to be of the order 
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of systematic errors expected from the use of isospin symmetry alone. We have looked for 
the origin of numerical differences. We studied both analytic forms of the distributions as 
well as dependence on numerical constants. 

Later, we have presented comparisons with further parameterizations, in particular 
the variants of TAUOLA RChL 2012, TAUOLA RChL and TAUOLA BaBar. All 
distributions were obtained with the help of Monte Carlo generations. We have prepared 
semi-analytical distributions to calculate such results as well. This is a necessary step 
toward htting with the help of multidimensional distributions, using variants of model 
independent methods presented in Ref. |1], but adopted to the case, when r leptons are 
relativistic (that is when momentum cannot be reconstructed easily). 

We have discussed limitations and constraints on use of models and symmetries like 
isospin symmetry resulting from insufficiently differential distributions of experimental 
data. We have pointed to large numerical consequences of interferences. We have also 
pointed out problems related to phase-space constraints for decay products of broad reso¬ 
nances. 

We should accept that predictions of the RChL models may agree with the data with 
precision not worse than ^ that is about 30% (in respect to total rate) only. Indeed, all 
of the presented parameterizations pass such condition and agree with each other at that 
level as well. Precision of the experimental data is at least one order of magnitude better. 
That is why, one has to concentrate the effort on arrangements for convenient comparisons 
of the data with model predictions. It is important that models are prepared in sufficiently 
flexible way (the necessary future adaptations can be introduced fast and correlations 
of parameters are understood) before the fits to experimental data are performed. Also 
experimental distributions need to be investigated, if they are detailed enough to constrain 
the underlying physics. 

We have discussed consequences of limited dimensionality of distributions used in data 
representations for results of fits. We could see in Ref. |T3], that the model of r —)■ 
decay modes, prepared only a year earlier |20] and fitted to two, one-dimensional, mass 
histograms, turned out to predict results for the invariant mas of pair rather poorly. 
It had to be updated. We are expecting this to repeat once three-dimensional experimental 
histograms are addressed. There is a factor of two discrepancy between RChL models and 
CLEO models results for some bins of the three-dimensional distributions (which were 
used by CLEO also as input to their fits). One-dimensional histograms may thus not 
be enough as experimental input for precise parameterization of hadronic currents. The 
three-dimensional ones of the type as used by CLEO [6] may be sufficient for r —)■ Sttz/,-. 
However, this conjecture may need to be checked before serious work on the r decay 
channels into 3 scalars with kaons, like that of Refs. ^51 06] is completed. Then, full 
complexity of hadronic currents is expected, and analysis similar to the one of Ref. [I] 
may be indispensable. Nonetheless, attempts to construct currents on the basis of fits 
to one-dimensional (three-dimensional CLEO style) distributions are of interest. It can 
demonstrate predictive power of models such as RChL beyond ^ precision level. On a 
practical side, once differences and options are understood in context of one-dimensional 
distributions, we are better prepared to fit the high precision data of today. 
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If operators like those in Ref. |1] are needed, then they have to be adopted for conditions 
of T leptons produced relativistically (making reconstruction of Uj. momentum non trivial). 
Our paper documents also a step of work into such direction, see also web page [47] . 

Even though the purpose of our paper is to review available options for hadronic cur¬ 
rents used in simulation of r decays into Svr, some speculations, why such choices were 
introduced in the past, were unavoidable. This is by far incomplete aspect of this work, 
and may be even misleading at some points. There is little documentation on internal 
discussions of experimental groups behind actual choice of analytical form of current pa- 
rameterizations used in data analysis and publications. Often complicated and sometimes 
unhnished activities, where questions of of great importance on systematic errors for mul¬ 
tidimensional distributions were risen, could not be referred. 

In general, we think experiments should be responsible for hts approving models quality. 
We hope our paper will contribute to discussion how relations between model developers, 
experimental physicists and people working on Monte Carlo and htting programs should 
look. 
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A Analytical distributions 

From the reference [13] documenting RCliL current and its £t to the experimental data it 
became obvious that having semi-analytical distributions available, simplifies and speeds 
up the htting procedure. As a part of this work we have prepared a program to calculate 
distributions from given hadronic current by integrating the structure function Wa (as 
dehned in Ref. [3|) within limits of a chosen bin. Program uses IntegralMultiple function 
of ROOT [38] libraries allowing acquisition of both one-dimensional and Dalitz distributions. 
We have checked the code to give proper results for TAUOLA CLEO, TAUOLA CLEO 
isospin intricate, TAUOLA BaBar, and TAUOLA RChL. Program offers easy way for 
re-fit of these models to the data of a form of three-dimensional distributions like in GLEG 
publication. Also, it is much easier, than if MG method is used, to perform analysis of 
statistical and systematical errors with the semi-analytical results free of statistical errors 
themselves. See Ref. [13] for details. 
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IntegralMultiple routine of ROOT library has several input parameters. In particular, 
minimum and maximum number of function evaluation requested for each bin. Those 
parameters affect precision and time of integration, therefore one has to hnd a balance 
between the two. We have found that integration has some stability problems if too low 
number of requested points is chosen. Taking into account time, precision, and stability of 
integration we choose minimum 5000 points and maximum 50 000 points per bin. With such 
setup it takes 6 sec (of CPU time on a typical PC processor) for TAUOLA CLEO and 16 
sec for TAUOLA RChL to get one-dimensional distributions and respectively, 47 and 120 
sec to get Dalitz plots with the same binning as presented in Appendix O One has to add 2 
sec for CLEO and 30 sec for RChL models for reinitialization whenever model parameters 
are changed. If higher precision is needed, doubling number of points per bin, requires 
50% more time for three-dimensional plots while it does not noticeably affect calculation 
of one-dimensional distributions. We expect that CPU time for Dalitz plots calculation can 
still be optimized. Fit with Dalitz plots is more desired, it gives better test of investigated 
model and has simpler systematics calculation; each individual event enters only one bin 
in contrary to one-dimensional histograms where it enters eacf@ of four histograms of 
invariant masses. In Ref. [H] optimization of htting procedure for TAUOLA RChL 
in case of one-dimensional histograms is presented. We hope that with similar approach 
we could prepare framework allowing fast htting to three-dimensional distributions for any 
theoretical model as well. Note, that current initialization of RChL is prepared for multiple 
channels of r decays, including not yet published ones with kaons |l5l Il6] . These, were not 
confronted with the data yet. For htting to 7r“7r“7r+ data, we can skip running routines 
needed for these other channels, to speed it up to the level which was used in Ref. |44] . 
We have not checked if IntegralMultiple will provide smooth distributions with model 
parameters variation. It is important for htting algorithms that derivatives with respect 
to model parameters are continuous, as we have found while working for Ref. jl4j . 

B Comparison of TAUOLA CLEO and Pythia CLEO 

In this appendix, we present comparison of alternative implementation of CLEO model 
for r —)■ as done in PYTHIA 8.201, with TAUOLA implementation. Thanks to 

comparison we can identify source of diherences and estimate its impact on generated 
event samples. 

By investigation of Pythia code, several diherences in parameterization were found and 
stored in Table El Parameters not mentioned in that table are the same as for TAUOLA. 
Hadronic current of Pythia CLEO is of the same form as described in Section 0] for 
TAUOLA CLEO isospin intricate. Internal arrangements of the code are slightly 
diherent. Due to this; in structure of the code some constants are folded with others. 
Parameter P^^ is one of such. It is not coded explicitly, but we can extract its value: 
0.784468371 GeV. 

Each event enters histograms for invariant mass of 7r“7r“, 7r“7r“7r''" and twice for 7r“7r“'' systems (as 
there are two possible pairs of 7r“7r+). 
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TAUOLA CLEO 

Pythia CLEO 

Mr = 1.777 GeV 

Ma, = 1.275 GeV 
= 0.700 GeV 
C3PI= 0.2384^ 

Mr = 1.77699 GeV 

Ma^ = 1.331 GeV 
is calculated, see section text 
C3PI = 0.238471.0252088 

/33 = BT3 = 0.37e-*‘^-^^^ 
(3^ = BT4 = 0.87e*o-537r 
(3^ = BT5 = 0.71e*°-®®’" 

/?3 = BT3 = 3.7 ■ 

13^ = BT4 = 8.7 ■ lO-7e*o-53^ 

(3^ = BT5 = 7.1 ■ lO-7e*o-56^ 


Table 6: Differences in numerical values of parameters used in TAUOLA CLEO and 
Pythia CLEO. 

Figures [H] and [15] of Appendix O represent Dalitz plots from TAUOLA and Pythia, 
respectively. Fig. [13 shows ratio of Fig. [T3to Fig. [TT] One-dimensional distributions of 
invariant mass of 7r~n~^, 7r“7r“, 7i~7;~7r~^ for both generators can be seen on Fig. [TD] There 
is systematic difference between generators. Disagreement may have appeared due to 
different form of hadronic current, or different numerical parameters. To test this, for both 
generators constants of Ref. [6] were taken. Agreement was not improved substantially, 
with the exception of invariant mass 7r“7r“7r+, where influence of ai mass and width turned 
out to be the main source of difference, as can be seen from Fig. [TT] In general, agreement 
only to about 10 % can be concluded except some less populated bins at ends of 
mass spectra. 




Figure 10: Invariant mass of TT tt’*', tt tt , tt n vr"''in TAUOLA CLEO (red) and Pythia 
CLEO (green). For details of histograms definition see caption of Fig. [T] 


It suggests that functional form of hadronic current is of higher numerical importance, 
than model constants. To test this, we have used TAUOLA CLEO isospin intricate 
current while having the same parameters in both generators. As in previous test, constants 
were set that of [6]. Resulting distributions of invariant masse^as well as their ratios can 
be seen in Fig. [12) nearly perfect agreement can be concludeco. 

To complete our investigation, comparison of TAUOLA CLEO isospin intricate 
and Pythia CLEO, with default initializations was performed. Generators use hadronic 
current of eqs. (|HD, 02]), and only numerical constants differ as those collected in Table EJ 

This is also true for Dalitz plots, see Fig. [T8| 
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Figure 11: Invariant mass of , 7T~7r~, in TAUOLA CLEO (red) and Pythia 

CLEO (green), when the same parameters are enforced for both generators. For details 
of histograms dehnition see caption of Fig. [H 




Figure 12: Invariant mass of in TAUOLA CLEO isospin intri¬ 

cate (red) and Pythia CLEO (green), when the same parameters are enforced for both 
generators. For details of histograms dehnition see caption of Fig. [1] 


Figure dn] shows resulting distributions of invariant masses and their ratios. One can see 
the biggest difference in invariant mass of due to distinct Oi mass and width, 

which are of highest importance among all numerical constants. 



Figure 13: Invariant mass of tt tt’*', tt tt , tt tt tt"'' in TAUOLA CLEO isospin intri¬ 
cate (red) and Pythia CLEO (green). For details of histograms dehnition see caption of 

Fig. [H 

Further, technical checks were performed: the same numerical constants and the same 
but simplihed hadronic current were installed in both generators. Those included: 

• reducing hadronic current to Fi of eq. ([5]) only (F 2 , F 3 = 0). 

• setting hadronic current to constant. 

• reducing hadronic current to Breit-Wigner of ai only. 
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In all these cases no discrepancies between results of TAUOLA and Pythia generators beyond 
statistical fluctuations were found. Samples of 3 million events were used. 

All considered, differences in produced by TAUOLA and Pythia distributions for the 
r —)■ 'K~'K~'K^i>r decay channel result only from parameters if the same variant of current is 
used in TAUOLA, otherwise generators perform identically. All differences are insignihcant 
in comparison to parameterization uncertainties. 

C Cleo style results of 3-dimensional distribution 

In this Appendix we collect 3 dimensional distributions for 3 scalar hnal states following 
dehnition used by CLEO [6]. For each hgure, eight Dalitz plots in si, S 2 invariant masses 
are given for restricted, respectively to 0.36-0.81, 0.81-1.0, 1.0-1.21, 1.21-1.44, 1.44-1.69, 
1.69-1.96, 1.96-2.25, 2.25-3.24 GeV^ ranges. 

Such representation does not constrain full differential nature of hadronic currents, but 
it was used in publication of experimental data of the most differential nature until now. 
It is an important reference point. 
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Figure 14: 8 Dalitz plots for slices in 0.36-0.81, 0.81-1.0, 1.0-1.21, 1.21-1.44, 1.44- 
1.69, 1.69-1.96, 1.96-2.25, 2.25-3.24 GeV^. Each Dalitz plot is distribution for TAUOLA 
CLEO in si, S 2 variables (GeV^ units). Si is taken to be the highest of the two possible 
values of , in each event. 

TT TT^ 
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Figure 15: 8 Dalitz plots for slices in Q^: 0.36-0.81, 0.81-1.0, 1.0-1.21, 1.21-1.44, 1.44-1.69, 
1.69-1.96, 1.96-2.25, 2.25-3.24 GeV^. Each Dalitz plot is distribution for Pythia CLEO 
in si, S 2 variables (GeV^ units). Si is taken to be the highest of the two possible values of 
I in each event. 

TT TT^ 


35 





































































































Figure 16: 8 Dalitz plots for slices in Q‘^: 0.36-0.81, 0.81-1.0, 1.0-1.21, 1.21-1.44, 1.44- 
1.69, 1.69-1.96, 1.96-2.25, 2.25-3.24 GeV^. Each Dalitz plot is distribution for TAUOLA 
CLEO isospin intricate in si, S 2 variables (GeV^ units). Si is taken to be the highest 
of the two possible values of , in each event. 

^ TT TT^ 
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Figure 17: Ratios of Pythia CLEO to TAUOLA CLEO Dalitz plots (Figs. [151IH!)- 
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Figure 18: Ratios of Pythia CLEO to TAUOLA CLEO isospin intricate Dalitz 
plots (Figs. [El ITHll . 
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Figure 19: Ratios of TAUOLA CLEO isospin intricate to TAUOLA CLEO Dalitz 
plots (Figs, [ini E]). 
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Figure 20: Ratios of TAUOLA BaBar to TAUOLA CLEO Dalitz plots in si, $2 
variables (GeV^ units). Consecutive plots correspond to slices in Q‘^\ 0.36-0.81, 0.81-1.0, 
1.0-1.21, 1.21-1.44, 1.44-1.69, 1.69-1.96, 1.96-2.25, 2.25-3.24 GeV^. 
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Figure 21: Ratios of TAUOLA RChL to TAUOLA CLEO Dalitz plots in si, S 2 variables 
(GeV^ units). Consecutive plots correspond to slices in Q^: 0.36-0.81, 0.81-1.0, 1.0-1.21, 
1.21-1.44, 1.44-1.69, 1.69-1.96, 1.96-2.25, 2.25-3.24 GeV^. 
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